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ABSTRACT 

In this paper the interaction problem of a through or a part- 
through crack with a stress-free boundary in a semi-infinite cylindrical 
shell is considered. It is assumed that the crack lies in a meridional 
plane which is a plane of symmetry with respect to the external loads 
as well as the geometry. The circular boundary of the semi -infinite 
cylinder is assumed to be stress-free. By using a transverse shear 
theory the problem is formulated in terms of a system of singular inte- 
gral equations. The line spring model is used to treat the part-through 
crack problem. In the case of a through crack the interaction between 
the perturbed stress fields due to the crack and the free boundary is 
quite strong and there is a considerable increase in the stress inten- 
sity factors caused by the interaction. On the other hand in the prob- 
lem of a surface crack the interaction appears to be much weaker and 
consequently the magnification in the stress intensity factors is 
much less significant. 

1. Introduction 

The main objective of this paper is to study the effect of a stress- 
free boundary on the stress intensity factors in a cylindrical shell 
which contains a part-through or a through crack near the boundary. From 
the corresponding elasticity solutions of flat plates it is known that 
when the stress field perturbed by a crack interacts with a free boundary 
there is a certain magnification in the related stress intensity fac- 
tors, implying that there would be a reduction in the resistance of the 
structure to fatigue crack propagation and fracture [1,2]. Experimental 
as well as analytical results also show that generally the stress 
intensity factors at that crack tip which is nearer to the free boundary 
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are greater than those at the other tip and consequently the crack would 
tend to propagate toward (and eventually perpendicular to) the free 
boundary [1-3]. This (nears that of the simple crack geometries in 
plates and shells the one which may be of the most practical interest 
would be the crack perpendicular to the free boundary. 

In shells additional magnification in the stress intensity factors 
is known to arise because of curvature [4-6], In this paper the com- 
bined effect of the stress-free boundary and curvature in a cylindrical 
shell is studied. It is assumed the plane of the free boundary is per- 
pendicular to the axis of the cylinder and the crack is in a meridional 
plane. Further, the external loads applied to the curved surfaces and 
the end of the cylinder are assumed to be such that the shear stresses 
a-j 2 and a ]3 vanish throughout the shell in the plane of the crack (Fig. 1). 
Consequently, in the crack problem the plane of the crack would be a 
plane of symmetry with respect to loading as well as geometry. 

2. Differential Equations and the Boundary Conditions 

The problem described in Fig. 1 is considered within the confines 
of linearized shallow shell assumptions and a Reissner typ*' ^nsverse 
shear theory [7]. Details of the formulation of a general snallow shell 
leading to a crack problem may be found in [4-6] and [8], The problem 
is conveniently formulated in terms of a stress function <f>, transverse 
displacement w and two auxiliary functions $ and a. Thus, referring 
to Appendix A for notation and for definitions of dimensionless and 
normalized quantities, the basic equilibrium equations of a cylindrical 
shell may be expressed as follows: 


(JL) 2 a n 

7 * ( T ; W 0 ’ 

(1) 

A, + a 2 a, 2 (1-kv 2 ) £$ = A-O-KV 2 ) $ , 

(2) 

<V 2 t l> - ip - w ~ 0 , 

(3) 
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where q(x,y) is tha transverse shear loading and t|> and a are related 

to rotations g and g by [5]: 
a y 

a - 3 jL + x( 1 "v) ft - llL k(1-v) 3£ / 1- 

" ax 2 ay ’ p y " ay 2 ax ' 

In terms of the stress function ^ the membrane resultants are 
defined as 


N = lit N - M a . 

xx W 2 - * yy ax 7 ^ * xy axay 


( 6 ) 


Similarly, the moment and the transverse shear resultants are given 
in terms of the displacement and rotations by 


M 


xx 


= fix* (_ 3T + v ly 4 * M yy = fif* (v 


ax ay 7 * 


M 


xy 


a Jbv f !?x . Zl) 

hT^ 2 'ay ' * 


ax 


(7) 


V = iw 

x ax 


V.. = 


_ aw 


ay 


+ 3, 


( 8 ) 


Before stating the boundary conditions one may note that the 
problem is linear and can be solved by using a superposition technique. 
Thus, the problem of a relatively thin-walled semi-infinite cylinder 
without the crack under given external loads may be considered separ- 
ately from the crack problem in which statically self-equilibrating 
crack surface tractions are the only external loads. In the former 
problem the external loads may be the transverse loading q and stress, 
moment, and transverse shear resultants applied to the end (x 2 =0) of 
the cylinder. In the plane of the crack x^=0 this solution is assumed 
to give only normal membrane and bending resultants and and no 
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shear stresses (i .e. , n 12 (o,x 2 ) 

a ^12^ ,X 2^ s y i ^ ,x 2^ a 

equal and opposite of N^-|(0,x 2 ) and M 11 (0,x 2 ) given by this solution 
are then used as the crack surface tractions to solve the crack problem. 
In this study it is assumed that the problem of the semi -infinite 
cylinder without the crack has been solved and N^(0,x 2 ) and M-|-j(0,x 2 ) 
are known functions. 

The part- through crack problem will be solved by using the line- 
spring model. The integral equations of this problem are obtained by 
modifying those of the through crack problem to account for the unbro- 
ken ligament. The main problem is, therefore, the through crack prob- 
lem which will first be formulated. Also note that because of symmetry 
it is sufficient to consider one half (x>0) of the shell only. Thus, 
referring to the insert in Fig. 2 and Appendix A for notation, equa- 
tions (l)-(8) must be solved under the following boundary conditions: 

Nyy(x»0) s 0, N^(x,0) = 0, M yy (x,0) - 0, 

(9) 

( 10 ) 
01 ) 
( 12 ) 

(13) 

(14) 


M xy (x,0) = 0, V y (x,0) = 0, (0<x<«) , 

N xy (0,y) * 0, M xy (0,y) - 0, V x (0,y) = 0, (-~<y<0) , 
N xx (+0,y) = F-j (y) , (-d^yob^ , 
u(0,y) * 0, (-»<y<-d 1 , -b 1 <y<0) , 

M xx (+0,y) = F 2 (y) , (-d^-bj) , 

3 x (0,y) = 0 , (-°°<y<-d 1 , -b^<y<0) , 


where F^ and F 2 are related to the uncracked shell solution by 

Nti(0,x ? ) ^^(OjXp) 

F l (y) = - . F 2 (y) a - > < x 2 =ay) - 


( 15 ) 
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3. Solution of Differential Equations 

In the perturbation problem under consideration q=0 and eliminating 
<P from (1) and (2) we obtain 


V 4 V 4 W + Aj 4 (1-kV 2 ) 



(16) 


Considering the symmetry of the problem, for the shallow cylindrical 
shell the solution of (IS) may be expressed as 


OO W 

w(x,y) - f T (x,a)e’ iya da + ~ f 2 (y, 3)cosx3 d@. (17) 

-00 Q 

Assuming the solution for f^ and f 2 of the form 


f^x.a) = R(a,m)e mx , f 2 (y,s) = S(e,n)e ny , 


(18) 


by substituting from (18) and (17) into (16) the characteristic equa 
tions for m and n may be obtained as 


m 8 -4a 2 m s +6a 4 m 4 -(4a 2 +KXj 4 )a 4 m 2 +a 4 [a 4 +X^ 4 (l+Ka 2 )J = 0 , (19) 

n 8 -(4$ 2 +xX 1 4 )n 6 +[63 4 +x 1 4 (l+ic8 2 )]n 4 -43 6 n 2 +3 8 = 0 . ( 20 ) 

The solutions of (19) and (20) have the property 

m j+4 = “ m j * Re ( m j) < °» (j“l*.*»»4) , (21) 

n j+4 * ' n j * Re ( n j) < °* • (22) 

Defining now 

R(a,nij) - Rj(oi) , S(3»nj) 15 Sj($) » (j = l»»**»8) , (23) 
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with (17) the solution of (16) satisfying the regularity conditions 
at x=« and y a -» may be expressed as 


4 m-x 

f» (x»a) = S R^(a)e J , (0<X«») , 

1 j«l J 

8 n,y 

f?(e»y) s i S . ( 0 )e J , (-«<y<0) . 
L j=5 J 


(24) 

(25) 


The functions R., (j=l,...,4) and S., (j s 5,...,8) are as yet unknown. 

J J 


Similarly, if we let 


4 >(x»y) = 2 ^r I g 1 (x,a)e’ ,ya da + - J g 2 (y,3)cosx3 d3 


-lya 


(26) 


noting that q=0 from (1) and (2) the functions g^ and g 2 may be obtained 
as follows: 


X-I 2 4 a 2 R.(a) m-X 

g^CXta) = -(y) S — y z — e J , Pj * mj 2 -ot 2 , (0<x<“) , 

J— 1 


j a 


j— 5 ^ 

Assuming the solution of (4) of the form 

oo oo 

-If 


a(x,y) = Tj^r h 1 (x,a)e“ iya da h £ h 2 (y, 3 )sinsx d3. 


(27) 


2 8 n . 2 S . (3) n.y 

g 2 (y»3) = (y) . s -^yr— ■ e J , s nj 2 -3 2 , (-~<y<Q). (28) 


(29) 


h-j (x,a) = A(a,r)e rx , n 2 (y,3) = B(s,s)e 


sy 


it can be shown that 

h-j(x,a) = A-j (a)e , r^ ~ -[a 2 + yyy] , (0<x<°°) 

h 2 (y»e) = B 2 ( 3 )e ^ , s 2 = +[3 2 + (-«<y<0) 

where A, (a) and Bjfs) are unknown functions. 


r l x 


(30) 


(31) 

(32) 
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The transverse shear theory used In this study Is equivalent to 
a tenth order system [4]-[6] (that Is, for example, the equilibrium 
equations of the shell can be expressed in terms of five second order 
differential equations in five "displacement" quantities u, v, w, $ 
and $ ) and the differential equations (1), (2) and (4) provide suffi- 
cient number of "integration constants" to account for all boundary 
conditions. For a unique solution of the problem, it is therefore 
necessary to use only the particular solution of (3) (after determining 
the function w). Thus, by assuming 


«(x,y) 


oo 


1 f 



-00 


0>l (x ,a)e ^ a da + ~ 


00 


0 


0 2 (y» 3 )cos 3 x d$, 


(33) 


the functions and satisfying the regularity conditions at x=« 
and y=-°° are found to be 


4 Rj(«) 

' j=l KP j ' 


m.x 
> j 


9 2 (y»3) « .2 


S^s) 


J-5 '‘V* 


, n j y 


(34) 


Ten unknown functions R., (j=l,...,4), S . , (j=5,..„,8), A, and 

J J I 

B 2 which appear in the solution given above are determined by using 
the eight homogeneous and two mixed boundary conditions (9 )- ( 1 4) . 


4. Derivation of Integral Equations 

By substituting now from the solution given in the previous sec- 
tion into (5)- (8) and by using the strain displacement relations 


c ij = I (u i,j +U j,1 +Z ,1 U 3.j +Z ,j U 3,l)' OJ * 1 > 2 > ’ 


(35) 


all membrane, moment and transverse shear resultants, the rotations, 
and in-plane displacements can be expressed in terms of the unknown 
functions Rj, S^, A-j and B 2 . These expressions are given in Appendix B. 
In (35) Z(xpX 2 ) gives the distance of the middle surface of the shell 
from the tangent plane. By substituting from Appendix B into (9) and 
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by Inverting the sine and cosine transforms the homogeneous end con- 
ditions may be expressed as follows: 


8 8 3 n, 2 i 

p — -4— s .(e) * J 


J-8 V ~ 3 


2ir 


4 AkRka) 
j=i d “ ’ 


(36) 


8 n. 3 . 

2 a 2 " ^<(8) 3 2n 

j=5 q j j 2n 


4 a 3 mjR.(a) 

I! / I ^ ^ t m. 1 ' *7 dot f 

j=i (m j +s ,p j 


(37) 


8 vB 2 -n ( 2 


2 i - rr) i K < 1 -' , > 2 r - n Ir) - Mr- '""-I 

jf 6 V e) + Mr~ s s 2 b 2 (s) - ' sr j >■,*, T^pMTO?+F] 


4 (vmj 2 -a 2 )m^R 1 *(a) 

1 kr -j 


i k( 1 - v) 2 ar^ 2 

+ -Jf^7+gzy- Aj(a)]da , 


(38) 


8 2611 , , r 4 21 aein f R < (a) 

jsg ^ V B) - ^z- 1 (s 2 2+ 3 2 )B 2 (e} ■ - ST kf, r^-TOjk 2 ) 


K(l-v)(a 2 +ri 2 )8 
+ 2(r 1 2 +e 2 ) A 7 (a) ]da 


(39) 


i< «F - -T- B ¥ B > - - 27 1% (^-ik^) 


(l-v)r, 2 / 

+ 2^7'Z+g2) A<| (a)]da . (40) 

Three more homogeneous conditions may be obtained from (10) which, with 

(36)-(40) may be used to eliminate eight of the ten unknowns. The 
mixed boundary conditions (11 )- (14) would then give a system of dual 
integral equations to determine the remaining two unknowns. However, 
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the problem may also be reduced directly to a system of integral equa- 
tions by introducing the following new unknown functions? 

G](y) a jy u(+0,y) , G 2 (y) « ~ 8 x (+0,y) , -»<y<0 . (41) 

From (12), (14) and (41) it may be seen that the conditions 
(12) and (14) are equivalent to 


Gj(y) B 0, -»<yodj , -b^<y<0, (j*l,Z) 

7 b i 

1 S,(y)dy = 0, (j=l,2) . 


(42) 

(43) 


We may now use (41) in place of (11)-(14) as two additional equations 
to express all unknowns (R^ , Sj, A-j and B 2 ) in terms of G-j and G 2 . Then 
(11) and (13) with (42) and (43) would provide the information to deter- 
mine G.j and G2. First, by substituting from equations (Bll) and (B9) 
of Appendix B into (41), it may be shown that 


- sr <r> : 


4 (l+v)a 2 -p. , 

Z jj—jr — Rj(a;a e y da , 

J - 1 J 


(44) 


s 2 (y) * - 


4 amjR^a) e _ iyo _ k(1.v)_ | „ 2fl y 0 

i- j=l V 


j da. 

(45) 


Then, by using the relations (B3), (B6) and (B7), after some manipula- 
tions it can be shown that equations (10), (44) and (45) are equivalent 
to 


4 m.R.(ot) 

2 -Vt-- 0 . 

j=l p j 


(46) 


4 m.R.(a) 1 

.1 -£ [S.~-T = 1 Wl-v)a + i]C 2 (a) , 


(47) 


j=l 
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(48) 


4 / V 

2 m.RAd) 3 
j*l J J 


i c 2 («) 


> 


4 

2 


j=l 


nijR^a) 


Ai 2 

{-i) C, (a) 


(49) 


A*j ( ot ) 3 2 Cg(a) , 


50) 


where 


Cj(a) 3 


-b. 


-d. 


Gjtyle 1 ^ dy , (j=l,2) . 


(51) 


form 


After solving (46) -(49) for and obtaining expressions of the 


-bj 

Rj(ct) 3 j* 2 Cj|^(ct»t)G^(t)dt » (j 3 l m * * *4) » 

a k=l 
“ d l 

from (50), (52) and (36)- (40) one can determine Sj and B 2 and find 
expressions of the form 
-b. 


(52) 


Sj (B) 3 


'1 2 


-d 

-b 


2 bj k (3,t)G k (t)dt , (j 3 5 . . ,8) 


1 


k=l 


1 2 


B 2 (6) 3 I 2 b 5k (e,t)G k (t)dt 

Id, ksl 


(53) 


(54) 


where Cj k and bj k are quite complicated but known functions. 

If we now substitute from (50), (52) -(54) , (Bl) and (B4) into 
(11) and (13) we obtain a system of integral equations which, with 
(43), would determine G.| and G 2> These equations are of the following 
form: 



•b, » 

1 2 


» 
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Tim ( 2 [f V H (x,a)e i ( t “ y ^ c< da + Y H (yyt,0)cQSX0d0]G.,(t)dt 

x-*+0 J. j«l l KJ Jo KJ 3 


B F k (y), (k«l,2), (~d-| <y<-b-, ) , 


(55) 


where V kj « and Y k ^ are known functions. The kernels defined by the 
functions V k j and Yy in (55) exhibit certain singular behavior which 
may be investigated by examining the asymptotic behavior of these 
functions for |ct|-*> and $-*». The process is relatively straightforward. 
First one needs to obtain the asymptotic expressions for m^, n . , r-j and 
s 2 defined respectively by (19), (20), (31) and (32) for large values 
of |a[ and 0 . These expressions may be obtained as follows: 



Then, by substituting from (56)-(59) into the expressions of Vy and 
Y kj ., the principal parts of the kernels can be separated and can be 
evaluated in closed form. After somewhat lengthy manipulations the 
integral equations (55) may be shown to reduce to 



( 1 1 , 6y 

^t-y t+y + (t+y)* 


rt^F )G i< t)dt 



2 

E 


j = 1 


k Tj . (t)dt = 27rF 1 (y), (-d 1 <y<-b 1 ) 


(60) 
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i k 2 .(y,t)G.(t)dt 
j«l d3 J 


2WA 1 * jF 2 (y), (-d 7 <y<-b 1 ). 


(61) 


where k^, (1,j*l,2) are known bounded functions and the loading func- 
tions F-j and Fg are given by (15). 

It should be pointed out that the dominant kernel found in this 
study » hich is given in (60) is -:entical to that found for the elas- 
ticity problem of a plane containing a crack perpendicular to a stress- 
free boundary [1], It should also be observed that the dominant 
kernels associated with the crack opening displacement and the crack 
surface rotation as given in (60) and (61), respectively, are identical, 
verifying the result found in [9] for flat plates. These are physi- 
cally expected results and mean that any peculiar behavior the shell 
solution may have for b=0 (i.e,, for an edge crack) must be identical 
to that of the cow^pandlng plane elasticity solution. 

For b>fi (fii> £ 2) the kernels of the integral equations (60) 
and (61) have simple Cauchy type singularity and consequently the 
solution is of the following form [10]: 

n k (y) 


G k (y) « 


[-(y+^Xy+b,)]* 


, (-d 1 <y<-b 1 ), (k-l ,2) 


(62) 


1 


where the unknown functions H-j and Mg are bounded in the closed interval 
[-dj, -b^], In this case the problem may easily be solved numerically 
by using a Gauss-Chebyshev type integration formula [11]. For b s 0 
the problem is an edge crack problem, the functions G 1 and Gg are 
bounded at. y=-b^=0, and equations (43) are no longer valid. In this 
problem too the Integral equations may again be solved by using a 
Gauss-Chebyshev integration formula and by, for example, assuming that 
H k (0) a 0 (in (62)) in place of conditions (43). 
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5. Stress Intensity Factors 

In the linearized shell theory used In this study the stresses 


are given in terms of the stress resultants as follows: 

CT Tj ( x 1 » Xg ) a Njj (Xi ,Xg) » ("I *2) , (63) 

b 12x 3 

^1 j ^ X 1 * x 2 ,x 3^ K “fjl j ( x i »Xg ) ’ (i»j ss 1»2) , (64) 

° r 3j( x l » x 2 * x 3 ) a ^j( x 1 ,x 2^”4t/^ ( J-l »2 ) , (65) 

where the "1n-p1ane ,( stresses may be combined as 

a ij^ x l ,x 2 * x 3^ = a ij (^i >Xg) ^ CT ij ^ X 1 ,x 2 ,x 3^ * 0 ><5*1*2) • (66) 


In the syrmietric problem under consideration is the only nonzero 
stress component in the plane of the crack x,=0. Thus, the crack 
problem is one of Mode I and ^(Xg^Xg), (x 2 ^-b,-d) is the only nonzero 
stress Intensity factor which, at the crack tips x 2 = -b and Xg=-d may 
be defined as follows: 


M-b.x,) * lim /2Tx7+6T 

J Xg-^-b c 


^11 ( 0 » Xg ,Xg ) , 


(67) 


k-j ( -d ,x^) - lim r/~2( Xg+d ) o , ^>|(0,Xg,Xg) • (68) 

Xg-v-d 

From (63), (64) and (66)-(68) it is seen that k-j would be a linear 
function in Xg and, analogous to (66), it would be convenient to express 
it as 


k^(Xg,Xg) « ^(xg) + jg kjjfxg), Xg = (-b,-d), - ^ < Xg < ^ , 


- 1 3 - 


(69) 


where and k b are associated with the "membrane" and "bending" 
stresses or and around the crack tips. Observing now that the 
left hand sides of the integral equations (60) and (61) give the expres- 
sions for N xx and M xx in the outside as well as the inside of the 
crack on the Xj a 0 plane and that only the terms with the dominant ker- 
nels would contribute to the singular behavior of N xx and M xx> the stress 
intensity factors k m and k b may easily be evaluated in terms of the 
unknown functions Gj and G 2 . Thus, by using the properties of singular 
Integrals [10], it may be shown that 

k m (-b) . - § ^ H, (-b, ) , k b (-b) = H 2 <-b,) , (70) 

k m ( -d ) = f & Hj(-dj) , kjjf-b) = £ /a ) . (71 ) 


6. The Part -Through Crack 

In this paper to solve the part-through crack problem shown in 
Fig. 1 the line spring model [12] is used. The details of the analysis 
with regard to the adoption of the model to Reissner plate or shell 
theory and the comparison of the results with the existing three- 
dimensional finite element solutions may be found in [13] and [14]. 
Following [14], the Mode I stress intensity factor along the crack 
frong is expressed as 

k-| (y ) = fTF [ff(y)g t (s) + m(y)g b (s)L ? = L(y)/h , (72) 

where the functions a and m represent the membrane and bending resultants 
along the net ligament -d^<y<-b-j and 9 t /S$ and (s=L/h) are the 

shape factors obtained from the plane elasticity solution of a strip 
having a thickness h, containing an edge crack of depth L and subjected 
to uniform tension or bending away from the crack region. Following 
are relatively accurate expressions of g t and g b : 
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g t U) = (1.1216 + 6.52005 2 - 12.38775 4 + 89,05545 s 

- 188,60805 s + 207. 38705 10 - 32.05245 12 ) » (73) 

g m (5) = /5 (1.1202 - 1.88725 + 18.01435 2 - 87.38515 s 

+ 241.91245 4 - 319.94025 s + 168.01055 s ) . (74) 


To obtain o(y) and m(y) first the integral equations (60) and (61) 
are modified as [13], [14] 

y y -bj 

-y t t (y)j g t (t)dt + r tb (y)j G 2 (t)dt + ^ j 


“di 


-d-i 




-d-,<y<-b 1 , 

Y bt (y)|^ G,(t)dt-Y bb (y)| s 2 (t)dt + j [(^ 


(75) 


l 


-d, 


-di 


‘ Wy + TZ+yF ‘ Tt+yP )G 2 (t >“ 


-b 


1 


2 k 2j .(y,t)G J .(t)]dt=F 2 (y) 
:d l j=1 


-d 1 <y<-b 1 (76) 

where the functions y- , (i,j s t,b) are given in [14] and the upper 
(i.e., -) and lower (i.e., +) signs are to be used for the outer and 
the inner surface crack, respectively. After solving (75) and (76) a 
and m are then obtained from 
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o(y) ■ E(y tt U i Y tb e x ) * m(y) = 6E(y bt U + Y bb e x ) , 


u(+0,y) 


(t)dt , 3 x (+0,y) - 



-d-j<y<-b^ » 


(77) 

(78) 


where + and - signs are to be used for the external and the internal 
surface crack, respectively. 


7. Numerical Results and Discussion 

For the through crack problem the calculated results are shown 
in Tables 1-4 and Figures 2-11. In the results given in Tables 1-4 
it is assumed that in the crack region the shell is subjected to uni- 
form membrane loading or bending, i.e., for the uncracked shell we 
have 

N-j •] ( 0 # Xg ) = Ng , *j (0 ® » (79) 

or 

N-| i ( 0 » Xg ) = 0 » M'j i (8»x 2 ) = Mq • (80) 

Thus the input functions of the integral equations are (see (15)) 

F v (y) = -N 0 /hE, F 2 (y) = 0 or Fj(y) = 0, F 2 (y) = -M Q /h 2 E . (81) 


In the tables the membrane and bending components of the stress inten- 
sity factor defined by (69) are normalized as follows: 


ml J a f /a ’ bl J o,^ 


, O. = (-d,«b) , i=(m,b) , (82) 


where refers to the crack tip and for uniform membrane and bending 
loads the corresponding stresses are given by 
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(83) 


«m 3 yn . % s 6 V h2 * 

The results are given for various combinations of R/h , a/h and c/a 
which are the characteristic dimensionless length parameters of the 
shell (see the insert in Fig. 2). The normalized stress intensity 
factors given in the tables for c/a * 1 are obtained from the edge 
crack solution. A partial explanation of the steep rise in the stress 
intensity factors at x 2 S! - d as c-*a may be found in [9]. For the membrane 
loading when R-*» or a->0 the bending component of the stress inten- 
sity factor tends to zero and k^ approach the corresponding flat 
plate solution given in [1]. The tables show the general magnifying 
effect of the decreasing crack distance c from the boundary and the 
curvature on the stress intensity factors. 

Some results obtained from the uniform external loads (79) and 
(80) are also shown in Figures 2-6. For a fixed relative crack dis- 
tance c/a = 1.5 Fig. 2 shows the effect of R/h and a/h on the membrane 
component of the stress intensity factor. As a/h->0 the results 
approach the corresponding flat plate values obtained in [1] which 
are also shown in the figure. Figures 3-6 show the edge crack results. 
Except for Fig. 5, qualitatively these results are quite similar to 
those obtained for the infinite shell problem given in [4]. Again it 
is interesting to observe that as a/h->0 k mm approaches 1.586 given 
by the plane elasticity solution for an edge crack regardless of the 
value of R/h. Even though this is the physically expected result, as 
in the limiting cases given in [9], the surprising aspect of it is 
its prediction by a shell theory. In an infinitely long cylinder 
having an axial line crack the curves corresponding to those given in 
Fig. 3 would tend to k^-1 as a/h->0. It is therefore tempting to look 
for a shape factor which, as in flat plates, would be independent of 
the crack length. However, multiplying the results given in [4] cor- 
responding to those in Fig. 3 for a/h>0 by 1.586, it may be shown that 
the relative magnification in k mm in the edge crack is considerably 
greater than that in the infinite cylinder. 
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It may be seen from Table 3 and Fig. 5 that as in the infinite 
cylinder case the bending stress intensity factors found for the 
edge crack too decreases with increasing crack length. In this problem 
the steep variation of k bb and its convergence to the plane elasticity 
solution for a/VO may seem to be somewhat unusual but should be physi- 
cally expected. In fact this behavior is nearly Identical to the 
bending solution of a flat plate with an edge crack given in [9]. 

In the Gauss-Chebyshev quadrature formulas used in solving the 
integral equations (60) and (61) the integrals on the left-hand sides 
are evaluated for certain discrete set of values y^ of normalized 
variable y. Therefore, one may easily obtain some type of Green's 
functions for the stress intensity factors by simply taking one element 
of the column matrices F-| (y^ ) and F 2 (y 1 >) non-zero at a time. This 
was routinely done and the results were printed in solving the problem. 
By using these results one can calculate the stress intensity factors 

k mm’ k bm* k mb* anct k bb any arb * tra Ty external loads 1^(0, x 2 ) and 
M-j i (0 »x 2 ) or F-j(y) and F 2 (y), For example, Figures 7-11 show some 
results for a concentrated wedge loading of a cylinder containing an 
axial edge crack. Figures 7 and 8 show the membrane and bending com- 
ponents of the Mode I stress intensity factor at the crack tip x 2 =-d 
k-|(-d,x 3 ) = kj(x 3 ) in a cylinder subjected to concentrated wedge forces 

N n(°’ x 2> a ‘ P6 ^ (84) 

shown in the figure. One may note that, as in wedge-loaded flat plate 
problems in the shell problem too the stress Intensity factor is pro- 
portional to l//a and as a/h->0 (or R-*») again the shell results con- 
verge to the plane elasticity values k^ = 0.58 (Fig, 7) and k bm = 0 
(Fig. 8). Figures 9 and 10 show the results analogous to Figures 7 
and 8 for membrane wedge forces applied at x 2 = -a, i.e,, for 

Ni i (0 » x 2 ) = -P<s(x 2 +a) . (85) 
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The effect of the wedge force location or the membrane component 
of the stress intensity factor for some selected shell and crack dimen- 
sions is shown in Fig. 11. The corresponding stress intensity factor 
in an infinite plate containing a semi -infinite crack and subjected 
to wedge forces P at a distance s from the crack tip is known to be 


k 


, P 
1 ’ n hvT 


( 86 ) 


Thus, since k^ becomes unbounded as s->0 in Fig. 11 P/h/s" rather than 
P/h/a is used as the normalizing stress intensity factor. It is seen 
that as s->0 the stress intensity factor in the shell approaches the 
flat plate value given by (86), As in the previous limiting cases 
discussed in this section in this problem too if the applied loads are 
very near the crack tips or the length of the (through) crack is very 
small compared to all other dimensions of the structure (including its 
thickness), then locally the stress state is indistinguishable from 
the corresponding plane elasticity solution. 

The stress intensity factors obtained for a part-through crack 
are given in Tables 5-12. Even though the crack profile L(y) can be 
any smooth function, the results are calculated only for a semi-elliptic 
inner or outer surface crack in the cylinder. That is, in solving the 
integral equations (75) and (76) it is assumed that (Fig. 1) 

I(x 2 ) = L 0 Ji-(~) 2 - 1 0 , 

y s x 2 /a, c-|=c/a, -d<x 2 <-b, -c^-l<y<-c^+l . (87) 


The problem is solved by assuming that the applied loads are given by 
(79) or (80) and hence the input functions in (75) and (76) are given 
by (81). In both cases only the Mode I stress intensity factor k-j 
is nonzero along the crack front. k t and k^ shown in the tables cor- 
respond to the normalized Mode I stress intensity factor obtained for 
uniform "tension" (79) or for "bending" (80), respectively, k^ is a 
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function of x 2 only (Fig, 1) and for the two loading conditions con- 
sidered is normalized as follows: 


k t (x) 


•<l(x 2 ) 


f * 


X ? +c L 

, X = — -y + c,. E 0 = f 


( 88 ) 


k b ( *> " 5FT 


k-j (x 2 ) 




x 2 +c L 

x = — =y + Cl , « 0 -f 


(89) 


where L Q is the maximum depth of the crack (Fig. 1) and the shape func- 
tions g t and g b are given by (73) and (74), One may note that the 
normalizing stress intensity factors used in (88) and (89) are those 
obtained from the plane strain solution of an infinite strip having a 
thickness h and containing an edge crack of depth L Q which is subjected 
to uniform tension N Q or bending M Q . For the crack depths L Q jsed 
in this study the shape factors g^ and g b are calculated to be 


5 0 * W" 

0.2 

0.4 

0.6 

0.8 

^ g t ( e 0 > 

1.3674 

2.1119 

4.035 

11.988 

9 b ( V 

1.0554 

1.2610 

1.915 

4.691 


Tables 5-9 show the normalized stress intensity factors k fc and k b 
for a certain range of R/h, a/h, c/a and L Q /h which are the dimension- 
less characteristic length parameters of the problem. It is assumed 
that the part-through crack may be on the outside or the inside sur- 
face of the cylinder. As in the through crack solution, in all cal- 
culations the Poisson's ratio v of the material is assumed to be 0.3. 
The previous studies varying v between 0 and 0.5 have shown that its 
effect on the stress intensity factors is relatively insignificant 
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(see, for example, [5]). Hence, the results given should be valid for 
most structural materials. The tables show that the stress intensity 
factor at the maximum penetration point of a semi -elliptic surface 
crack appears to be relatively insensitive to the crack distance c from 
the free end. 

Some sample results showing the variation of the stress intensity 
factor along the crack front are given in Tables 10-12. It is seen 
that generally at a given crack depth the stress intensity factor at 
the location nearer to the stress-free boundary is greater than that 
at the corresponding point which is farther. However, this skewness 
or non-symmetry in the distribution of the stress intensity factor with 
respect to the midpoint x = 0 or ^ - -c of the crack seems to be 
surprisingly small. 
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APPENDIX a original page 19 

OF POOR QUALITY 

Dimensionless and Normalized Quantities 

x * X * y 8 T » 2 * T » b l * ? » c i a £ . d ! = | . 

(A.l) 

u » u^/a , v * u 2 /a , w * u 3 /a , 

(A. 2) 

. F(X-| »X« ) 

3 X 3] , 3 y 3 2 , tf)(x,y) - a 2 Eh , 

(A. 3) 

a xx s a ll /E » a yy = a 22 /E » a xy s ~1T * CT xz * "IT * ^yz = ' 

3 > (A. 4) 

Nl 7 N«« N-.« m N. • 

N xx = Tif ’ N yy " TiT ’ N xy * "RT > °ij = "R 1 > < 1=1 > 2 ) • 

(A. 5) 

m M 11 ^22 ^19 h 6 M 4 j 

M xx * m ’ M yy * fiZf - M xy = fPf • a ij = Tfz 1 • (U-1.2) 

(A.6) 

v x V hB » V y “ V hB » 8 " T2TT+v) 

(A. 7) 

4 * 12(l-v 2 )aVh 2 R 2 , = 12(l-\> 2 )a 2 /h 2 , K . E/BX 4 , 

(A. 8) 


For dimensions and notation see Figure 1. 
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ORIGINAL PAGE 19 
OF POOR QUALITY 

APPENDIX B 

Expressions for stress resultants, displacements and rotations. 

r w 4 r” Q 

N xx (x,y) * z c^Tjjdct + j | n j 2 LjjCOS8xd3» 


N yy (x,y) * - 


-00 


N xy (x,y) = -1 


Z a 2 fTlj 2 Tjjda - 


E a 3 mjT-jjda + 


CO 

f 8 

Z L 1: ji3<*cos0xde , 
0 5 


00 ^ 

| 0 njL-|jSinexd 8 , 

0 


oo 

r 4 r 8 

M xx (x,y) * j s(mj 2 -va 2 )T 2 jda - E( 3 2 *vnj 2 )L 2 jCOS 3 xd 0 


(l-v)l 


a riT 0 da + 




1 3 aa T TiF 


8 s^L^cosexda 


00 
t 4 


00 

r 8 


Myy(x,y) = s(vmj 2 -a 2 )T 2 jda - E(v8 2 -nj 2 )L 2 jC0S8xd8 


- ^ f « r lT ,da - ^ 


J M '13 
-00 


00 
f 4 


Si 0 COS 0 Xd 0 , 


00 
f 8 


M^(x,y) = -i (l-v)j s a mjT 2j .da - (l-v)J E 3 n j L 2 j sinsxd3 


J 

-00 


w no 

- f (a^r^jTqda + (s 1 2 +3 2 )L 3 sin3xd0 , 


1 ' 3 


mx’ j v3 i 

0 


00 QQ 

V x (x,y)=^( \ PjtnjTj.jda - ^ { | MjLySlnwi* 


00 00 

- 1j 0 Tjda + [ s^jSlnsxdS . 


( 8 . 1 ) 

(B.2) 

(8-3) 

(8.4) 

(B.5) 

( 8 . 6 ) 

(B.7) 
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UHIGINAL PAGE 13 

op poor quality: 


V y (x,y) 


1 Khx 4 


TO 

( f° 


P J T 2j da * Kj f 


to 

r o 

5 q j n j4j COS0XdS 

0 


« co 

- | I^Tjda - j 8 LjCOSSxdS , 


B x (x,y) - ^ 


» ^ 09 

s m j T 2j da - J | 0 Lgjsinexde 

-oo o 


- i I a Tgda + 


TO 

1“ 

-TO 

IhX 4 


TO 


s-| singxdg , 


3y(x»y) s - a 


oo a> 

f 4 hx^ f 8 

2 a Tjjda + — j S HjLjjCOSSxdS 

-TO 0 


o 

■i 


r l T 3 d “ ‘ 


g Bg cosgxdg 


CO TO 

u(x,y) = | 2 [(2+v)a 2 - m j 2 ^ m j T ij do1 + | | (2+v)gL 1 jSingxdg 


+ (y) 2 x w(x,y) , 


TO TO 

v(x,y) = -if z(va 2 +m. 2 )aT ? .da - £(g 2 +vn , 2 )L, . dg , 

li j cj J s J i J n.* 

-oo 0 J 


where 


1 A, 2 Ri (a) m.x-iay 

Tn(X.y,Ci) - 4- T7--4-7- e i , 


id 1 **'’ 0 ' - ^rir-pTz- 


, R.-(a) m.x-iay 

T 2j ( x.y.a) = Jprc 


(B.8) 

(B.9) 

(B.10) 

(B.n) 

(B.12) 
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Tg(x >y I®} 


L^fy.s) * 

L 2 j(y.B) * 


, , •’I*"' 1 ''*' 

1 ii JrW“ )e 


§ £S£*w v . 


n,y 


2 3 1 p / „ \ ^ j 

?HxT^T s j (s)e 


sii^i B 2 (e)e S2y 


( y * 3 ) ~ 


Table 1. Membrane component of the normalized stress Intensity factors 
In a cylindrical shell containing an axial through crack 
along -d<X 2 <-b. The end of the cylinder (x2*0) Is stress- 
free and the crack surfaces are subjected to uniform membrane 
loading (Nji»»-N 0 , Mii»0). 
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10 6.092 4.654 
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Table 2. Bending component of the normalized stress intensity factors 

in a cylindrical shell containing an axial through crack along 
-d<X 2 <”b. The end of the cylinder X2 IB! 0 is stress-free and 
the crack surfaces are subjected to uniform membrane loading 
(Nji— N 0 , M n -0). 
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Table 3. Bending component of the normalized stress Intensity factors 

in a cylindrical shell containing an axial through crack along 
-d<X2<-b. The end of the cylinder x 2 "0 is stress-free and 
the crack surfaces are subjected to uniform bending (N^»0, 
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Table 4. Membrane component of the normalized stress Intensity factors in a 
cylindrical shell containing an axial through crack along -d<X 2 <-b. 

The end of the cylinder X2*0 is stress-free and the crack surfaces are 
subjected to uniform bending (Njj»0, Mj)“-M ). 
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0.151 

0.171 

0.169 


10 

-0.511 

-0.824 

-0.771 

-0.840 

-3.223 

- 1.190 

-0.847 

-0.826 

-O .838 


1 

0.092 

0.053 

0.053 

0.072 

0.024 

0.034 

0.043 

0.046 

0.071 

10 

2 

0.195 

0.126 

0.133 

0.133 

0.021 

0.087 

0.126 

0.126 

0.135 

3 

0.237 

0.197 

0. 186 

0.174 

-0.057 

0.131 

0.160 

0.165 

0.175 


10 

0.068 

-0.166 

-0. 140 

-0.192 

-1.886 

- 0.369 

-0.206 

-0.160 

-0.191 


2 

0.118 

0.070 

0.073 

0.087 

0.036 

0,046 

0.060 

0.071 

0.087 

25 

3 

0.173 

0.120 

0.126 

0.131 

0.034 

0.080 

0.106 

0.115 

0.130 


10 

0.312 

0.210 

0.164 

0.145 

-0.769 

0.079 

0.118 

0.155 

0.145 


1 

0.036 

0.018 

0.014 

0.013 

0.058 

0.020 

0.014 

0.011 

0.011 

100 

2 

0.056 

0.029 

0.028 

0.039 

0.024 

0.023 

0.023 

0.026 

0.037 


3 

0.085 

0.046 

0.046 

0.060 

0.026 

0.034 

0.039 

0.051 

0.060 


10 

0.254 

0.224 

0.205 

0.196 

0.033 

0.154 

0.185 

0.189 

0.198 


1 

0.035 

0.016 

0.004 

0.018 

0.062 

0.020 

0.014 

0.005 

0.015 

200 

2 

0.039 

0.020 

0.016 

0.023 

0.040 

0.017 

0.014 

0.014 

0.021 

3 

0.056 

0.030 

0.029 

0.042 

0.025 



0.025 

0.041 


10 



0.158 

0.163 

0.042 


0.138 

0.146 

0.164 
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Table 5. Normalized stress intensity factors in a cylindrical 

shell with a stress-free end containing an axial semi- 
elliptic surface crack and subjected to uniform local 
membrane loading or bending, R/h ■ 5. 


a 

£ 

L - 

0 

0,2h 

L o- 

0 . 4h 

L * 

0 

0.6h 

L * 
0 

0.8h 

h 

a 

k t (0) 

W 

k^foT 

W 

k £ (0) 

W 

W 

W 




Outer Crack, R/h « 5, 

V * 0. 

3 



l.l 

.843 

.833 

.561 

.507 

.292 

.189 

.090 

.008 

1 

1.5 

.835 

.824 

.543 

.485 

.275 

. 1 69 

.083 

.000 


2 

.832 

.820 

.537 

.477 

.2/0 

.162 

.082 

-.001 


10 

.827 

.815 

.529 

.468 

.266 

.157 

.081 

-.002 


1.1 

.904 

". 898 

.687 

.651 

74To 

.329 

.142 

7 062“ 

2 

K5 

.900 

.894 

.678 

.639 

.400 

.315 

.137 

.056 


2 

.898 

.891 

.673 

.633 

• 395 

.309 

.135 

.054 


10 

.894 

.887 

.664 

.623 

.387 

.299 

.133 

.051 


1.1 

.930 

.926 

.757 

.730 

.497 

.429 

.191 

.116 

7 

1.5 

.927 

.922 

.748 

.719 

.486 

.415 

.185 

.108 

J 

2 

.925 

.920 

.744 

.714 

.480 

.408 

.182 

.104 


10 

.922 

.917 

.736 

.705 

.472 

.397 

.179 

.100 


1.1 

. 966 

7954 

75Ti 

7555 

7551 

.638 " 

.342 

.275 

10 

1.5 

.964 

.962 

.866 

.849 

.675 

.627 

.334 

.265 

2 

.963 

.961 

.863 

.846 

.670 

.620 

.330 

.260 


10 

.961 

.959 

.857 

.838 

. 660 

.607 

.322 

.249 





nner Crack, R/h “ 5, 

v * 0.' 

) 



1.1 

T83I 

7520 

.539 

.480 

.274 

7168 

.085 

.001 

1 

1.5 

.821 

.809 

.517 

.454 

.257 

.145 

.079 

-.006 


2 

.816 

.803 

.508 

.443 

.250 

.137 

• 077 

-.008 


10 

.809 

.795 

.497 

.429 

.244 

.129 

.076 

-.010 


1.1 

7885 

.877 

7541 

.596 

.362 

7159 

.123 

.038 

2 

1.5 

.879 

.871 

.626 

• 578 

.347 

.250 

.117 

.031 


2 

.876 

.868 

.620 

.571 

• 342 

.243 

.116 

.029 


10 

.871 

.862 

.609 

.558 

.334 

.233 

.113 

.026 


l.l 

.911 

.905 

.703 

7557 

.428 

• 347 

7155 

.073 

7 

1.5 

.907 

.901 

.691 

.654 

.415 

.330 

.151 

.066 

J 

2 

.905 

.898 

.685 

.646 

.408 

.321 

.147 

.062 


10 

.901 

.894 

.676 

.635 

.399 

.310 

.145 

.058 


l.l 

.960 

.957 

.846 

7525 

7527 

.577 

7155— 

.215 

10 

1.5 

.958 

.955 

.839 

.820 

.616 

.562 

.277 

.203 

2 

.957 

.954 

.835 

.815 

.610 

.554 

.272 

.197 


10 

.955 

.951 

.829 

.807 

.600 

.541 

.266 

.189 
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Table 6. Normalized stress intensity factors In a cylindrical 

shell with a stress-free end containing an axial semi - 
elliptic surface crack and subjected to uniform local 
membrane loading or bending, R/h « 10. 


1 

a 

c 

L - 0.2h 
o 

L q ■ 0.4h 

L - 0.6h 
o 

L » 0.8h 
o 

h 


laaiHaa 

BIBBBS81 

WQ1IMQ1 

BHaiBSBi 


Outer Crack, R/h “ 10, v - 0.3 



T7T 

7551 

7 $30 

.555 

.^00 

.285 

7752 


. oo4~ 

1 

1.5 

.832 

.821 

.536 

.476 

.268 

.160 

ii 

-.003 

I 

2 

.828 

.816 

.528 

.467 

.262 

.152 

Hji a 

-.005 


10 

.823 

.811 

.521 

.458 

.257 

.146 

RE H 

-.007 


HI 

.900 

7894“ " 

7575 


.394 

.310 

BEH 

.050 

9 

Da 

.896 

.890 

.664 

.624 

.380 

.291 


.041 


D 

.894 

.887 

. 660 

.619 

.376 


■ 

.039 


turn 

.891 

.884 

.653 

.610 

.369 

.278 

.121 

.037 


TTT 

7923 

.923 


.718“ 


.405 

.172 

.095 

9 

1.5 

• 925 

.920 

Hr i \ 

.709 

Hpi 

.392 

.167 

.088 

J 

2 

.924 

.919 

m ' 

.704 

ffiSK 

.385 

.164 

.085 


10 

.921 

.916 

R 31 

.695 

.452 

.376 

.161 

.081 


1.1 

.970 

TsfcB 

m m 


mnam 

.661 

.351 

.289 

10 

1.5 

.968 

.966 

DU 



.651 

.344 

.280 

2 

.968 

.965 



vm* 

.646 

.340 

.275 


10 

.966 

.964 

MM 

.856 

HP* 

>1 Dggpjg 

.635 

.333 

.266 


Inner Crack, R/h * 10, v » 0.3 



in 

fin 


.542 

7555 

.275 

7759 

7655 

.000 

1 

Da 

EH 


.521 

.458 

.258 

.147 

.078 

-.006 

1 

mm 

mm 


.511 

.446 

.250 

.137 

.076 

-.009 


CM 

IBS 

.798 

.500 

.433 

.243 

.128 

.074 

-.011 


1.1 

W:f: :m 

R:l« 

. 646 

Km 

.363 

.271 

sn 


9 

1.5 


.874 

.629 

.582 

.345 

.248 


n 


2 

EES 

.871 

.622 

.574 

.339 

.240 

.no 



10 

mm 

.866 

.612 

.562 

.331 

.230 

.108 



in 

b m 

.908 

7765 

.669 

.423 

K» 

aa 


■a 

Da 

RSI 

.903 

.692 

.654 

.408 

Kill 

HP® 

.056 

J 

mm 

.907 

.901 

.687 

.648 

.403 

Dr 1H 


.053 


tm 

.903 

.897 

.678 

.637 

.394 

.304 

helM 

.050 


1.1 

79^3 

.961 

7553" 

.837 

.635 

7533 

.287 

.218 

1 n 

1.5 

.961 

.959 

.847 

.830 

.625 

.574 

.277 

.206 

1 U 

2 

.960 

.957 

.844 

.826 

.619 

.567 

.273 

.201 


10 

.958 

.956 

.839 

.819 

.611 

.556 

.267 

.194 


J 

A 


: * 
1* 


Y 
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Table 7. Normalized stress intensity factors in a cylindrical 
shell with a stress-free end containing an axial semi 
elliptic surface crack and subjected to uniform local 
membrane loading or bending, R/h ■ 25. 


a 

£ 

L * 0.2h 
o 

kb 

| 


L » 0 . 8h 
° 

h 

a 

W0BBIQB 

mm 

msai 

Laaiiaai 

BKHIBSSI 


Outer Crack, R/h - 25, v * 0.3 



T7T 

ttot- 

785“ 

.5$t ' 

“.w- 

"281 

7T7S 

7o57“ 

.001 

t 

1.5 

.829 

.818 

.530 

.470 

.263 

.154 


-.005 

i 

2 

.825 

.813 

.522 

.459 

.257 

.146 

.076 

-.008 


10 

.820 

.808 

.514 

.449 

.251 

.139 

.075 

-.010 

... 

l.l 

7888 

“S91 


.629 

T3S2 

.296 

.124 

.042 

0 

1.5 

.892 

.885 

.653 

.610 

.365 

.273 

.116 

.032 


2 

.890 

.883 

.647 

.603 

.359 

.265 

.113 

.029 


10 

.887 

.880 

.641 

.596 

.354 

.259 

.112 

.027 


1.1 

.925 

.920 

.735 

.705 

.455 

.381 

.157 



1.5 

.922 

.917 

.725 

.693 

.441 

• 363 

.149 

Rfj jff 

i 

2 

.920 

.915 


.688 

.436 

.357 

.147 



10 

.918 

.913 

• 715 

.681 

.430 

.350 

.145 

■RS « 


1.1 

.973 

.971 

7591 

758o 

T706 

.672 

.344 

.286 

10 

1.5 

.972 

.970 

.887 

.875 

.698 

.662 

.335 

.275 

2 

.971 

.969 

.885 

.873 

.695 

.657 

.332 

.271 


10 

.970 

.968 

.881 

.868 

.688 

.648 

.326 

.264 


Inner Crack, R/h - 25, v = 0.3 



ITT 

.835 

7825 

.543 

-MS— 

.275 ^ 

00 

vD 

7583 


1 

1.5 

.826 

.814 

.523 

.461 

.258 

.148 

E B 

-.007 


2 

.821 

.808 

.514 

.450 

.251 

.139 

.075 

-.010 


10 

.814 

.801 

.503 

.436 

.244 

.129 


-.012 


1.1 

.892 

7885 

.652 ' 

.611 

.367 

.276 

.119 

.034 

0 

1.5 

.885 

.878 

.635 

.590 

.348 

.252 

.110 

.024 

4 

2 

.882 

.874 

.627 

.580 

• 339 

.242 

.107 

.021 


10 

.878 

.869 

.618 

.569 

.332 

.231 

.105 

.018 


1.1 

• 917 

.911 

.711 

7677 

wmm i 

.345 

.144 

.061 

2 

1.5 

.912 

.906 


.660 

.408 


.135 

.051 

j 

2 

.910 

.904 

.691 

.653 

.401 

.314 

.133 

.047 


10 

ivrri 

.900 

.683 

.643 

.393 

.303 

.130 

.044 


1.1 

.965 

.963 

.859 

.844 

TTp 

.592 

727? 

.212 

1 n 

1.5 

.964 

.962 

.853 

.837 

.628 

.579 

.270 

.201 

1 u 

2 

.963 

.961 

.850 

.834 

.622 

.573 

.266 

.195 


10 

.961 

.959 

.846 

.828 

.615 

.563 

.261 

.189 
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Table 8. Normalized stress intensity factors in a cylindrical 

shell with a stress-free end containing an axial semi- 
elliptic surface crack and subjected to uniform local 
membrane loading or bending, R/h = 100. 


a 

c_ 

EB 

3.2h 

mvi 


mm 



h 

a 

sn 

BS&I 

mm 

ram 

W2M 

B1SI 

am 

Vor 




Outer Crack, R/h * 100, v * 

0.3 




mm 


mtm 

.547 

.490 

.276 

mmm 

7oF2 

-.002 

1 

D m 

® 

Wttm 

.528 

.467 

.260 

Bm m 

SiJriMi 

-.007 


mm 



.519 

.456 

.254 

IKK 

.075 

-.009 


19 

.818 

KM 

.509 

.444 

.247 

Wfim 


-.012 


mm 

.896 

m m 

7552 

7521 

.374 

.286 

.120 

.036 

0 

Da 

.890 

mm 

.646 

.602 

.356 

.263 

.112 

.027 


D 

.887 


.639 

.594 

.349 

.254 

.109 

.023 


EH 

.884 

.877 

.632 

.586 

.342 

.245 

.106 

.020 


i .i 

.923 

sum 

.726 

.695 

msssm 

.366 

HH 

.068 

•i 

1.5 

.918 


m is 

.680 

wTM 

.343 


.057 

J 

2 

• 917 

.911 

n Still 

.673 

.417 

.334 

aim 

.052 


10 

.915 

.909 

n 

. 667 

.412 

.327 

Wfim 

.050 


l.l 

.973 

.971 


.873 

.680 

.643 

.302 

.241' 

10 

1.5 

.972 

.970 


.869 

.673 

.635 

.296 

.233 


2 

• 971 

.970 

Epp 

.867 

.670 

.630 

.292 

.229 


10 

.970 

.969 

EB 

.863 

.664 

.623 

.288 

.223 


Inner Crack, R/h ** 100, v * 0.3 



1.1 

KHJB 

.825 

wEam 

.485 


. 1 66 

7o8T" 

nmm 

1 

1.5 

.827 

.815 

.525 

.463 

.258 

.148 

.077 

-.008 


2 

.822 

.810 

.516 

.453 

.252 

.140 

.075 

-.010 


10 

.816 

.803 

.506 

.440 

.245 

.131 

.073 

-.012 


1.1 

7594 

7 W“ 

7555 

.615 ” 

.368 

.278 

.117 

.033 

2 

1.5 

.888 

.880 

.640 

.595 

.351 

.256 

.110 

.024 


2 

.885 

.877 

.632 

.586 

.343 

.246 

.107 

.021 


10 

.881 

.873 

.623 

.575 

.334 

.235 

.104 

.017 


1.1 

.920 

.915 

.718 

7686 ""” 

.432 

.353 

.145 

~7o5r" 


1.5 

• 915 

.910 

.704 

.669 

.413 

.330 

.135 

.051 


2 

.913 

.907 

.697 

.661 

.405 

.319 

.131 

.046 


10 

.910 

.904 

.690 

.651 

.397 

.308 

.128 

.042 


l.l 

• 967 

.965 

7852 

".848 

7534 

T5W~ 

. 263 

“7194“ 

10 

1.5 

.966 

.964 

.857 

.841 

.624 

.575 

.254 

.183 

2 

.965 

.963 

.854 

.838 

.619 

.569 

.251 

.179 


10 

.964 

.962 

.850 

.833 

.612 

.560 

.246 

.173 
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Table 9. Normalized stress Intensity factors In a cylindrical 

shell with a stress-free end containing an axial semf- 
elllptlc surface crack and subjected to uniform local 
membrane loading or bending, R/h * 200. 








































Table 10. Distribution of normalized stress intensity factors along the 
crack front In a cylindrical shell with a stress-free end 
containing an axial semi-el 1 iptic surface crack and subjected 
to uniform local membrane loading or bending, x * (x 2 +c)/a, 
R/h » 10, a/h * 1, c/a « 1.1. 


BB 

1 




IS 

mm 


mm 

■89H 

mm 

M 

mm 

wm 


Outer Crack, R/h - 10, a/h ■ 1, c/a « 1 . 1 


.929 

.564 

— TS7A — 


" .615 


"Tpo 

.101 

.127 

.828 

.680 

.755 

.511 

.601 

.290 

.335 

.099 

.102 

.688 

.748 

.799 

.533 


.291 

.292 

.099 

.072 

.516 

.796 

.821 

.548 

.545 

.291 

.248 

.095 

.044 

.319 

.827 

.830 

.556 

.521 

.290 

.210 

.090 

mSSM 

.108 

.840 

.832 

.558 

.505 

.288 

.187 

.088 


0 

.841 

.830 

.555 

.500 

.285 

.182 

.087 

.004 

-.108 

.837 

.828 

.551 

.497 

.283 

.181 


.005 

-.319 

.816 

.818 

.537 

.498 

.277 

.192 

.085 

.017 

-.516 

• 777 

.800 

.516 

.505 

.267 

.217 

.086 

.034 

-.688 

.721 

.769 

.487 

.516 

.255 

.246 

.085 

.056 

-.828 

.648 

• 719 

.452 

.528 

.243 

.274 

.081 


-.929 

.553 

.639 

.414 

.533 

.232 

.306 

.078 

.097 


Inner Crack, R/h * 10, a/h * i, c/a =1.1 


.929 

.584 

.674 

HE 1 

.610 

.281 

• 371 

1119* 

.122 

.828 

.679 

.754 

1 

.594 

.282 

.325 

m* iifli 

.097 

.688 

.745 

.796 


.564 

.282 

.281 


.068 

.516 

.792 

.816 


.533 

.282 

.236 

.091 

.040 

.319 

.821 

.824 

.544 

.507 

.281 

.198 

.087 

.019 

.108 

.834 

.824 

.544 

.489 

.277 

.174 

.085 

.003 

0 

.834 

.823 

.542 

.484 

.275 

.169 

.084 


-.108 

.830 

.820 

.538 

.481 

.273 

.168 

.083 


-.319 

.809 

.811 

.524 

.483 

.266 

.179 

.082 

.043 

-.516 

.772 

.794 

.504 

.491 

.257 

.204 

.082 

.029 

-.688 

.717 

.765 

.477 

.504 

.246 

.234 

.081 

.051 

-.828 

.646 

.717 

.444 

.519 

.235 

.264 

.077 

mSm 

-.929 

.551 

.638 

.408 

.526 

.225 

.297 

.075 

HU 
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Table 11. Distribution of normalized stress intensity factors along the 
crack f**ont in a cylindrical shell with a stress-free end con- 
taining an axial semi-elliptic surface crack and subjected to 
uniform local membrane loading or bending, x ■ (xg+cl/a. 


L o /h 

0.2 

o.4 

0.6 

o.8 

* 

X 

im 

m 

mm 

n 

■BHk 

— 

SHHH 

k b 



Outer Crack, R/h 

m 10, a/h * 3, c/a * 1.1 



.929 

.549 

.639 

.456 


.330 

.457 

.151 

.203 

.828 

.671 

.749 

.541 

.654 

.372 

,458 

. 161 

.109 

.688 

.768 

.824 

.616 

.683 

.411 

.454 

,170 

.169 

.516 

.845 

.876 

,678 

.703 

.444 

.440 

.174 

.142 

.319 

.898 

.908 

.723 

.7 15 

.460 

.421 

.173 

• 117 

. 108 

.925 

.922 

.745 

.719 

.475 

.408 

.172 

.098 

0 

.928 

.923 

.746 

.713 

.475 

.405 

.172 

.095 

-.108 

.924 

.921 

.742 

.716 

.472 

.404 

.171 

.096 

-.319 

.894 

.903 

• 713 

.704 

.455 

.408 

.169 

.111 

-.516 

.838 

.868 

.662 

.683 

.425 

.416 

. 165 

.130 

-.688 

.759 

.814 

.593 

.654 

.384 

.419 

.156 

.149 

-.828 

.659 

.736 

.512 

.617 

.337 

.412 

.141 

.163 

-.929 


.625 

.423 

.563 

.289 

.403 

.126 

.171 


Inner Crack, R/h ■ 10, a/h ■ 3 ( c/a ■ 1.1 


.929 

.554 

Mzsm 

Mj ~ 

.616 


.459 

"iur 

.195 

.828 

.675 

.753 

.546 

.659 


.450 

.150 

.176 

.688 

.768 

.824 

.610 

.675 

.392 

.430 

.154 

.148 

.516 

.840 

.870 

.659 

.680 

.411 

.400 

.154 

.116 

.319 

.888 

.896 

.691 

.678 

.422 

.368 

.150 

.089 

.108 

.911 

.907 

.704 

.672 

.425 

.346 

.148 

. 068 

0 

.913 

.908 

.704 

.669 

.423 

.341 

.148 

.064 

-.108 

.909 

.905 

.699 

.666 

.419 

.340 

.146 

.065 

-.319 

.881 

.889 

.675 

.659 

.406 

.347 

.143 

.080 

-.516 

.829 

.858 

.633 

.648 

.384 

.365 

.141 

.100 

-.688 

.754 

.808 

• 574 

,631 

• 354 

.380 

.135 

.122 

-.828 

.658 

.734 

.503 

.606 

.318 

.387 

.124 

.142 

-.929 

.537 

.626 

.421 

.560 

.278 

.388 

.114 

.155 
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Table 12, Distribution of normalized stress Intensity factors along 
the crack front In a cylindrical shell with a stress-free 
end containing an axial seml-el 1 iptlc surface crack and 
subjected to uniform local membrane loading or bending, 
R/h * 25, a/h « 1, c/a * 1.1, x * (x 2 +c)/a. 


V h 

0.2 

0.4 

0.6 

0.8 


■9H 

n 

■Si 

ia 

mm 

mm 

■m 

k b 

X 

Outer Crack, R/h 

■ 25 , a/h * l , c/a * 1 . 



.929 

.581* 

.675 

.479 

.614 

.285 

.377 

.099 

.125 

,828 

.680 

.755 

.509 

.599 

.286 

.331 

.097 

.099 

.688 

.747 

.798 

.531 

.570 

.287 

.287 

.096 

.070 

.516 

.795 

.820 

.545 

.541 

.287 

.242 

,092 

.04 1 

.319 

.825 

.829 

.553 

.517 

.286 

.204 

.087 

.019 

.108 

.839 

.830 

.553 

.500 

.283 

.181 

.085 

.004 

0 

.839 

.829 

.551 

.495 

.281 

.176 

.084 

.001 

-.108 

.835 

.826 

.547 

.492 

.278 

.176 

-.084 

.003 

-.319 

.814 

.816 

.534 

.494 

.272 

.187 

.083 

.014 

-.516 

.776 

.799 

.513 

.502 

.263 

.212 

.083 

.032 

-.688 


.769 

.485 

.513 

.252 

.242 

.084 

.054 

-.828 

WiUm 

.719 

.450 

.526 

.240 

.271 

.083 

.076 

-.929 


.639 

,413 

.532 

.230 

.304 

.079 

.096 


Inner Crack, R/h 

* 25, a/h * 1 , c/a * 1 . 



.929 

.583 

.674 

.475 

.609 

.280 

.369 

.096 

.121 

.828 

.678 

.753 

.504 

.593 

.281 

.324 

.094 

.096 

.688 

.745 

.796 

.525 

.563 

.281 

.280 

.094 

.066 

.516 

.792 

.817 

.539 

.534 

.281 

.235 

.090 

.039 

.319 

.822 

.825 

.545 

.508 

.280 

.196 

.086 

.017 

.108 

.835 

.826 

.546 

.491 

.277 

.173 

.084 

.002 

0 

.835 

.824 

.543 

,486 

.275 

.168 

.083 

-.001 

-.108 

.831 

.822 

.539 

.483 

.272 

.168 

.082 

.000 

-.319 

.811 

.812 

.526 

,485 

.266 

.179 

.081 

.012 

-.516 

.773 

.795 

.506 

.494 

.258 

.205 

.082 

,029 

-.688 

• 719 

.766 

.479 

.507 

.247 

.235 

.081 

.051 

-.828 

.646 

.718 

.446 

.521 

.235 

.265 

.077 

.074 

-.929 

.552 

.638 

.410 

.528 

.226 


.075 

.093 
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Figures 


Fig. I 
Fig. 2 

Fig. 3 

Fig. 4 

Fig. 5 
Fig. 6 
Fig. 7 

Fig. 8 

Fig. 9 
Fig. 10 
Fig. 11 


Geometry of the semi-infinite cylinder with a part-through 
crack. 

Variation of the membrane component of the normalized stress 
Intensity factors In a seml-lnflnlte cylinder containing a 
through crack and subjected to uniform membrane loading 
N n -N 0 In the crack region. k mrn (-b)*l .204 and k mm (-d)**l .097 
shown In the figure are the stress Intensity factors In a semi- 
Inftnlte plate with' the same crack size, relative crack location 
and loading as the shell. 

Membrane component of the normalized stress intensity factor 
In a semi - inf Ini te cylinder containing an axial edge crack and 
subjected to uniform membrane loading NjpNq In the crack 
region. 

Bending component of the normalized stress intensity factor in 
a semi -Infinite cylinder containing an axial edge crack and 
subjected to uniform membrane loading Njj^Nq in the crack 
region. 

Bending component of the normalized stress Intensity factor in 
a semi-infinite cylinder with an axial edge crack which is sub- 
jected to uniform bending Mjj*M 0 in the crack region. 

Membrane component of the normalized stress intensity factor in 
a semi-infinite cylinder having an axial edge crack and subjected 
to uniform bending Mjj*M 0 in the crack region. 

Membrane component of the normalized stress intensity factor 
in a semi-infinite cylinder containing an axial edge crack and 
subjected to concentrated membrane wedge forces P at the cir- 
cular boundary. 

Bending component of the normalized stress intensity factor 
In a semi-infinite cylinder having an axial edge crack and sub- 
jected to concentrated membrane wedge forces P at the circular 
boundary. 

Membrane component of the normalized stress intensity factor 
in a semi-infinite cylinder having an axial edge crack and sub- 
jected to membrane wedge forces P atthe midpoint of the crack. 

Bending component of the normalized stress intensity factor in a 
semi-infinite cylinder having an axial edge crack and subjected 
to membrane wedge forces P at the midpoint of the crack. 

Effect of the location of concentrated membrane wedge forces on the 
membrane component of the normalized stress intensity factor in 
a semi-infinite cylinder containing an axial edge crack. 
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,204 

,087 



Variation of the membrane component of the normalized ‘.tress 
intensity factors in a semi- inf in i te cylinder containing a 
through crack and subjected to uniform membrane loading 
N n =N 0 in the crack region. k mrT 1 (-b) = l .204 and l< mni (- i d) ! =i f 097 
snown in the figure are the stress intensity factors in a semi- 
infinite plate with the same crack size, relative crack location 
and loading as the shell. 
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Fig. 3 Membrane component of the normalized stress intensity factor 
in a semi-infinite cylinder containing an axial edge crack and 
subjected to uniform membrane loading N|]=N 0 in the crack 
region. 
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Bending component of the normalized stress intensity factor in 
a semi-infinite cylinder containing an axial edge crack and 
subjected to uniform membrane loading N^=N 0 in the crack 
gio 
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Fig. 5 Bending component of the normalized stress intensity factor in 
a semi-infinite cylinder with an axial edge crack which is sub 
jected to uniform bending M|j=M 0 in the crack region. 
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Fig. 6 Membrane component of the normalized stress intensity factor in 

a semi - i nf i n i te cylinder having an axial edge crack and subjected 
to uniform bending Mjj=M 0 in the crack region. 


SEM 


Membrane component of the normalized stress intensity factor 
in a semi-infinite cylinder containing an axial edge crack and 
subjected to concentrated membrane wedge forces P at the cir- 
cular boundary. 





F the normalized stress intensity facto 
cylinder having an axial edge crack and 
ted membrane wedge forces P at the circi 
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Membrane component of the normalized stress intensity factor 
• n ? semi-mf inite cylinder having an axial edge crack and sub- 
jected to membrane wedge forces P atthe midpoint of the crack. 
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Effect of the location of concentrated membrane wedge forces on the 
membrane component of the normalized stress intensity factor in 
a semi-infinite cylinder containing an axial edge crack. 



